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Abstract: 

Hybrid metal-organic perovskite solar cells have astounded the solar cell community with 
their rapid rise in efficiency over the past three years. Despite this success, the basic pro¬ 
cesses governing the photogeneration of free charges, particularly their dynamics and effi¬ 
ciency, remain unknown. Here we use ultrabroadband pulses of THz frequency light to see 
the intrinsic photophysical properties of single crystal lead halide perovskite just femtosec¬ 
onds after a photon is first absorbed. Our spectra reveal the dynamics and efficiencies of 
free charge creation, the remarkable ease in which they move through the lattice and the 
complicated interplay between free and bound charges in these materials. 

The stellar rise of hybrid metal halide perovskites displaying impressive photovoltaic and 
water photolysis efficiencies combined with a composition of earth abundant elements, has 
lead to a flurry of spectroscopic work probing their underlying photophysics. Only very recently 
have large single crystals become available, allowing a glimpse into their intrinsic optical and 
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electronic properties A key spectral region for these materials is 1 - 100 meV, encompass¬ 
ing proposed exciton binding energies that determine charge photogeneration efficienc)^®, carrier 
scattering rates governing charge transport, and lattice vibrations whose role in screening the exci¬ 
ton remains controversial jn this work, we apply time-resolved multi-terahertz spectroscopy 
to probe photoexcitations over this entire spectral range in single crystal CHsNHsPbla with sub- 
100 fs resolution. We observe ultrafast dynamic screening of lattice phonons by free charges with 
remarkably high charge mobilities of 500-800 cm^/Vs. Optical pumping at the band edge, where 
exciton dissociation occurs purely through coupling to the thermal bath, reveals an exciton bind¬ 
ing energy of 49 ± 3 meV. At high photon flux, a dynamic exciton Mott transition is observed as 
mobile charges screen the insulating exciton gas leading to dissociation into a conducting plasma 
on a sub-1 ps time scale. Our measurements pave the way for a new understanding of fundamental 
excitations in this material, providing a direction toward optimization of this promising class of 
solar cells. 

The overall power conversion efficiency of organo-metal halide perovskite solar cells has 
climbed to 20.1% in a span of only three years the fastest rate of increase compared to all 
other solar cell technologies. This fast pace in device development has been matched by an effort 
to understand the intrinsic charge generation, recombination and transport properties of these ma¬ 
terials. Among the key questions remaining are the exciton binding energy the exciton/free 
charge generation dynamics the energetics of lattice excitations and their role in exciton disso¬ 
ciation Time-resolved terahertz spectroscopy (TRTS) is a powerful technique to probe these 
low energy excitations immediately after photoexcitation on femtosecond time scales Previous 
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THz studies on the trihalide and mixed halide perovskites have foeused on frequeney-averaged 
THz absorption, probing energies below 10 meV^^E^E^. These measurements on thin film samples 
have plaeed initial lower limits on the mobility, shown suppressed bi-moleeular reeombination, 
and indieated the role of the exeiton binding energy in the generation of free eharges on a ~ 2 
ps time seale Here we perform eomplete energy and time-resolved ultra-broadband THz spee- 
troseopy in the 8-100 meV speetral range with 40 fs temporal resolution after photoexeitation on a 
single erystal of the hybrid metal halide perovskite CHsNHsPbIs. Photoexeitation at the band edge 
probes the kineties of mobile eharge generation via dissoeiation of exeitons without the influenee 
of injeeted exeess kinetie energy. Our data reveal the generation dynamies and quantum yields of 
free eharges via thermal dissoeiation governed by an exeiton binding energy of 49 ± 3 meV. We 
further observe remarkable earrier mobilities that reaeh as high as 800 em^/Vs on ultra-short time 
seales, and the dynamie sereening of the exeiton at high photon flux. 

Fig 1(a) shows the sehematie of the optieal pump - THz probe measurement of the large, 
mm^ sized single erystal of CHsNHsPbIs, probing the (100) faeet of the erystal with image pro¬ 
vided. The unit eell of the eubie phase is also shown, with the average position of the CH3NH3 
eation indieated by the orange ball in the eenter of the unit eell. In these experiments, an opti¬ 
eal pulse photoexeites the sample and the indueed eharge degrees of freedom are then probed by 
transmission or refleetion of a phase-stable, single eyele eleetromagnetie transient with frequeney 
eomponents in the THz range, ealled a THz pulse. The ineident THz pulse (Fine), is shown in Fig. 
1(b) along with the transmitted (Etrans) and refleeted (Erefl) pulses. The ultra-broadband THz pulse 
Einc(t) eontains speetral eomponents from 1-30 THz (4-125 meV or 33-1000 om~^), as shown 
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Figure 1: Multi-THz spectroscopy of CHsNHsPbla (a) A schematic of the time-resolved 
multi-THz spectroscopy experiment is given, where a 40 fs, 1.56 eV pump pulse photoexcites 
the CHsNHsPbIs crystal at normal incidence, colinearly with a multi-THz probe pulse with 
Fourier spectral components spanning 1-30 THz (4 - 125 meV). The sample is a single crys¬ 
tal CHsNHsPbla approximately 1.5 x 1.5 mm in size with the (100) facet as the probed surface, 
(b) In the absence of a photoexcitation, the incident (black), transmitted (blue) and reflected (red) 
multi-THz pulses after interaction with the crystal are shown with (c) corresponding Fourier am¬ 
plitude spectra. The amplitude of the pump-induced differential reflected field, \AE{u, 600/s) |, is 
also shown in green, (d) The dark state dielectric function probing the (100) plane in transmission, 
with areas of no information indicated by the hatched regions. 
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in the corresponding Fourier amplitude spectrum in Fig 1(c). While the reflected pulse retains all 
incident Fourier components with noticeable dispersion, Etrans(t) is strongly chirped with no spec¬ 
tral components below 25 meV due to absorption and reflection arising from an infrared active 
phonon and dielectric relaxation of the CH 3 NH 3 cationic. The change of amplitude and phase of 
the THz pulse after transmission provides the complex dielectric function i(u;) = €i(cl>) -f * 62 ( 0 ;) 
via inversion of the Fresnel equations, and is shown in Fig. 1(d). The spectra is dominated by an 
Pb-I bending mode at ~ 8 meV previously observed in Raman and THz spectroscopy GUCSl jn ad¬ 
dition, a broad GHz dielectric contribution has also been identified from microwave measurements 
due to the relaxation of the CH 3 NHj|' cation 1^. The CH 3 NH 3 group is free to rotate and undergo 
librations in the room temperature tetragonal phase giving a Debye contribution to the dielectric 
function with a picosecond relaxation time. As a result of these two lattice responses, the opti¬ 
cal transmission is strongly suppressed for energies below 25 meV. The real dielectric function, 
ei(co), responsible for renormalizing the Coulomb interaction of injected electrons and holes and 
therefore the exciton binding energy Eb is nearly constant at 5.5 for 40 - 100 meV range. 
Below 40 meV near resonant with the exciton binding energy, however, ei (ce) is quite dispersive. 
The dynamics of exciton-phonon coupling should therefore be considered when discussing the 
appropriate dielectric function^. 

We now turn to the time-resolved THz spectra after photoexcitation. The reference multi- 
THz waveform in the absence of optical pumping, Eref(t), is shown in Eig. 1(b) and Eig. 2(a) 
after normal incidence reflection from the ( 100 ) crystal facet, showing slight reshaping due to 
the dispersion of the crystal. The power spectrum of this pulse, shown in Eig. 1(c), contains 
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frequency components spanning the entire 30 THz bandwidth. The pump-induced change in the 


reflected THz electric field AE(t,rp) =Epump(f, Tp)-E^ef(t) at a given pump-probe delay time Tp is 
measured at variable delays producing the two-dimensional time-domain data set shown in Eig. 
2(b). The Eourier amplitude of the modulated field in quasi-equilibrium is given in Eig. 1(c), 
showing spectral information is obtained between 8-80 meV. The significant phase shift compared 
to the reference pulse indicates a strong change in the dispersion after charges are photo-injected 
into the crystal. This is a signature of injected mobile charges screening the lattice from the applied 
field, within the pump pulse penetration depth of 1 /xm in the crystal The amplitude of the 
AE(t,rp) response signifies the onset of photoconductivity in the sample, and reaches its maximum 
within 400 fs after photoexcitation, seen in the constant Tp waveforms in Eig. 2(c). We note that 
this is 10 times the temporal resolution of our technique and thus mobile carriers are not generated 
instantaneously but rather through exciton dissociation^. 

The corresponding spectral amplitude and phase of the pump-induced change in THz re¬ 
flectance, defined as Af{u, Tp) = AE{u, Tp)/Erefiuj), are shown in Eig. 3(a) and (b), respectively 
for pump-probe times from 0 - 400 fs at a pump fluence of 80 /xJ/cm^. The data can be understood 
by the analytic formula^ 



( 1 ) 


where a = + ia 2 is the complex ac photoconductivity, h = n + in, is the complex crystal index 

of refraction and tq is the complex crystal reflection coefficient given by fo = {1 — n)/{I+ n). In 
the absence of photoconductivity, Af{uj) = 0 as well as in the limit for a highly reflective crystal 


with n » 1 where tq —1. Eig. 3(c) shows the complex index of refraction of the crystal used 
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Figure 2: THz transients and pump-induced differential waveforms (a) The refleeted THz 
transient from the unexeited CHaNHaPbla erystal (100) faeet and (b) the two-dimensional AE(t,rp) 
= Epump(t,rp)-Eref(t) data set for pump-probe delay times Tp < 1 ps. (e) Seleeted euts in the map 
along eonstant Tp are shown, representing the buildup of photoeonduetivity in the erystal. 
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to describe fo, composed of a picosecond Debye dielectric relaxation and Lorentzian response 
modelling the 8 meV infrared active phonon in agreement with Fig. 1(d). The high energy index 
of 2.35 is chosen to agree with our transmission data of Fig. 1(b) (Coo ~ 2.35^ = 5.52) and low 
energy limit from literature ~ 30® The corresponding reflectivity amplitude |ro| is shown in 
Fig. 3(c) exhibiting a sharp drop to a minimum just below 20 meV. Thus the screening of the lattice 
by injected mobile charge carriers is responsible for the large Ar response observed in Fig. 3(a), 
and by Kramers-Kronig the corresponding phase change in Fig. 3(b). By inversion of Eq. 1, the 
complete a{u, Tp) spectra can be extracted in temporal slices following photoexcitation, revealing 
the underlying ultrafast charge transport characteristics. 

An example of the a{u,Tp) spectrum recorded at = 1 ps after excitation is shown in 
Fig. 3(d) by inversion of Eqn. 1. The data is very well described by a simple Drude model 
a{oj) = 1 -Lt ’ "'here is the mobile charge density and r is the momentum scattering time 
determining the charge carrier mobility through the relation /i = er/m* with m* = 0.2me being 
the charge carrier effective mass^^. The Drude response confirms the band nature of transport in the 
single crystal CHaNHaPbIs as opposed to charges in a disordered film where hopping conductivity 
occurs. The scattering time r can be estimated simply by the crossing point of cxi and <72, occurring 
at an energy Ec = h/(27rr) ~ 10 meV. This provides an estimate of r ~70 fs corresponding to 
a remarkably high /i ~ 620 cm^/Vs for a solution processed semiconductor. This measurement 
represents the highest mobility observed in these materials to date although recent improvements 
in grain size and large crystal fabrication have seen /r on the order of 150 cm^/Vs.^ Moreover, the 
amplitude of the photoconductivity directly yields the mobile carrier density ~ 3.3 x 10^^ cm“^. 
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estimated by the Drude de limit adc = noe^. The mobile eharge earrier generation effieieney is 
given by no/nph ~ 12% where Uph = 2.9 x 10^® em“^ is the absorbed photon density given 
the pump fluenee of 80 /rJ/em^ and the 10% refleetion of the pump at 795 nm^^. This result is 
eonsistent with the expeeted result of the Saha-Langmuir equation^ where a mobile earrier yield 
of 13% is expeeted given a 50 meV exeiton binding energy at 293 K and an exeiton effeetive mass 
niex = 0.11 me^. 

The results of global fitting of the Drude model to the eomplete two-dimensional time- 
frequeney ^{u) maps, fitting both real and imaginary eomponents simultaneously, are shown in 
Fig. 4. The dynamies of the mobile eharge generation effieieney for two pump fluenees is shown 
in Fig. 4(a), exhibiting a slight transient enhaneement for Tp < 10 ps at the lower fluenee. The 
Saha-Langmuir equation aeeounts for this behaviour predieting a deerease in effieieney to 10% for 
the elevated fluenee due to the greater probability of eleetron-hole bindingl^. The remaining ~ 90% 
of exeitations, however are bound exeitons that are unable to sereen the lattiee direetly due to their 
eharge neutrality. If pump photons exeite earners above the band edge, however, exeiton dissoeia- 
tion will be infiueneed by the initial exeess kinetie energy, providing higher mobile eharge yields. 
The goal of this study, however, is to probe the energeties of thermal dissoeiation governed solely 
by the exeiton binding energy and so we leave this to future work. The mobility is also extraeted 
through the Drude fits as /i = er/m* and is found to be relatively fluenee independent, as shown in 
Fig. 4(b). A 20% deerease is observed in earrier mobility on a sub-pieoseeond time seale, whieh 
we eonjeeture eould be due to the eharging of trap states or some eharge-mediated phonon seatter- 
ing meehanism. Despite the slight reduetion in the transient mobility, the steady-state mobility of 
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Figure 3: Free charge carrier screening of IR phonons and ultrafast Drude conduction, (a) 

Amplitude and (b) phase of the differential multi-THz refleetanee shifted by steps of 0.5 for elarity. 
(e) The model of the eomplex index of refraetion of the unexeited erystal used to deseribe the 
substrate dispersion and extraet the eomplex photoeonduetivity. (d) The eomplex eonduetivity 
speetrum with Drude speetrum fits reeorded at 1 ps following exeitation. 
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Figure 4: Free charge generation, charge carrier mohilities and exciton dissociation, (a) The 

mobile charge carrier generation efficiency defined as the ratio of the charge density to the ab¬ 
sorbed photon density, at pump fluences of 80 and 140 /rJ/cm^. (b) Charge carrier mobility derived 
from the 80 and 140 /iJ/cm^ data sets with global complex fits to the Drude model, (c) The two- 
dimensional data set of the differential reflectance amplitude at a fluence of 810 /iJ/cm^, showing 
the early time region where a low energy transient peak appears at 12 meV. Vertical lines indi¬ 
cate cut regions shown in (d) at 12 meV and 20 meV corresponding to exciton and free charges, 
respectively. The dynamics of exciton dissociation are represented by the decay of the 12 meV 
peak concomitant with the 1 ps free charge density risetime. 
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550 cm^/Vs is orders of magnitude higher than other solution proeessed semiconductors. 

The induced mobile charges, however, should also screen the exciton as discussed in a recent 
observation of band filling effects and so one expects more efficient generation at very high 
carrier densities. The critical Mott density that defines the threshold for the efficiency of this 
screening, however, has been estimated at Ucru ~ 1 x 10^® cm“^, nearly an order of magnitude 
higher than densities in Fig. 3. To investigate the effects of charge screening of the exciton, a two- 
dimensional data set was taken at F= 810/rJ/cm^ with the |Ar(a;)| map shown in Fig. 4(c). The 
response is not well described by a Drude model alone, as it exhibits a transient low energy peak at 
12 meV, below the main reflection peak at 20 meV. The dynamics of this excitation are shown by 
the cut in the map along the Tp axis for the two peaks. The onset of the lattice screening occurs on a 
time scale of 0.5 - 1 ps, larger than the 40 fs temporal resolution of the measurement and indicating 
the charge generation is driven by exciton dissociation. The 12 meV peak initially tracks the rise 
time of the free charge 20 meV peak, however decays to the noise level by the time the lattice 
response is fully screened at Tp = 1 ps. This indicates a coupling between the two excitations and 
we conclude that the 12 meV peak is due to an intra-excitonic orbital transition given the binding 
energies for trap states have been recently estimated to be much larger, on the order of 100 - 400 
meV from the band edge.^^We note if excitation into trap states was the dominant photoexcitation 
at the band edge, one would not expect a Drude conductivity spectrum but rather a conductivity 
which is activated at a characteristic frequency. Even if the quantum efficiency reaches 4.5% 
predicted by the Saha-Langmuir equation at this elevated fluence, ud ^ 1.3 x 10^® cm“^ which 
exceeds Ucrit- In this regime, one expects all excitons are dissociated forming an uncorrelated 
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plasma. The disappearance of the 12 meV peak is therefore due to the dynamic screening of the 
injected exciton population as the free charge density increases beyond the critical Mott density on 
a sub-picosecond time scale. 

In conclusion, we have performed optical pump - multi-THz spectroscopy in the energeti¬ 
cally relevant 8 - 100 meV spectral region on a single crystal lead halide perovskite. Our spectra 
reveal the dynamic screening of infrared active phonons and background Debye relaxation of the 
CH 3 NH 3 cation, due to the photogenerated free charge population. Sub-picosecond conductivity 
spectra directly reveals free charge densities and remarkably high mobilities of 800 cm^/Vs on 
sub-picosecond time scales. Charge generation proceeds via exciton dissociation which is com¬ 
plete after 1 ps, with quantum yields of free carriers governed by a 50 meV binding energy. A 12 
meV signature of the exciton is observed, dynamically screened on a time scale comparable to the 
generation of mobile carriers and at a carrier density consistent with the excitonic Mott transition. 

Methods 

Sample preparation Single-crystals of CH 3 NH 3 Pbl 3 suitable for characterization were grown 
from a solvent mixture comprising of aqueous HI (57% w/w, 5.1 mL) and aqueous H 3 PO 2 (50% 
w/w, 1.7 mL). In a typical procedure, a 20 mL scintillation vial was charged with the colorless 
solvent mixture and heated to boiling (ca. 120° C). Addition of solid PbO (670 mg, 3 mmol) and 
CH 3 NH 3 CI (202 mg, 3 mmol) led initially to the formation of a black precipitate which rapidly dis¬ 
solved leading to a clear bright yellow solution. The boiling hot solution was then capped with sev¬ 
eral layers of parafilm tape to ensure that the composition of the vapors was maintained throughout 
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the crystal growth process. Omission of this step leads to the crystallization of CH3NH3Pbl3.H20, 
instead. On standing, upon reaching ambient temperature a countable number (in the range of 
10 -100) of small CH 3 NH 3 Pbl 3 crystals begin to form, which act as initial seeds for the subse¬ 
quent crystal growth from the supersaturated supernatant solution. Well-formed, faceted crystals 
of rhombic dodecahedral crystal habit were obtained after 2 weeks with crystal sizes ranging from 
1-4 mm among several batches, with the size of the crystals being inversely proportional to the 
number of the initial crystal seeds. The crystals were collected manually by decanting the mother 
liquor, pressed dry with a soft filtration paper and thoroughly dried under a stream of N 2 gas. 


Characterization The 795 nm, 40 fs pulse output of a Ti:sapphire regenerative amplifier is split 
into three beams. Charge carriers are injected the pump beam resonant with the band edge at 
1.56 eV to minimize initial excess energy, thereby suppressing direct exciton dissociation through 
this channel The penetration depth of the pump pulse is taken to be 1 /rm determined by thin 
film absorption, and it is assumed this is the same for a single crystal. Phase stable THz pulses 
are generated via a two-colour laser plasma in dry air and the electric field is directly detected 
after reflection off the crystal facet by an air-biased coherent detection scheme with 30 THz band¬ 
width 1^. Both Ej.e/(t) and AE(t,rp) are measured simultaneously through a double modulation 
scheme, required to minimize the influence of systematic timing variations. All measurements 
were performed under dry air environment and no evidence of oxidation was observed during the 
experiments. The background complex refractive index of the crystal h{uj) = and the 

dielectric function e{uj) is given by the Eorentzian phonon and Debye components 


eiu) — Coo + 




+ 


(cUq — — iuj'f 1 — iujTo 


( 2 ) 
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The high frequency dielectric function Coo = 5.5 determined from transmission measurements in 
Fig. 1(b) and the Pb-I bending mode a;o/27r = 1.96 THz and the Debye relaxation time td = 
2.5 ps. The spectra were found to be relatively insensitive to the phonon relaxation rate, set to 
a reasonable 7 = 0.1 THz. The Saha-Langmuir equation, used to determine the exciton binding 


energy from the x = riD/nphot ratio, is defined as^^ 
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